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Abstract 

In the framework of Quark-Gluon-String Model developed recently in ITEP we 
calculate spectra of charmed particles D, D s , A c , H c , £l c in hadron-hadron collisions 
taking into account the decays of 5-wave resonances like D* , D*, S c , £*, H*, E' c , 
and Q*. We describe the bulk of the existing data on D, D* , and A c production 
in irp and pp collisions and predict the yield of charmed particles in Ti~p and 
reactions at hyperon beam energies of 340 GeV/c and 600 GeV/c. Because of significant 
production of baryon resonances our predictions for unfavored fragmentation differ 
from predictions of other models which do not take resonance production into account. 
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1 Introduction 



The physics of charmed particles, e.g. particles which occupy the place in between 
light hadrons (hadrons containing u—, d- and s-quarks) and heavy hadrons (with 6-quark), 
is a subject of a great interest. On one hand, the methods available for light hadron sector: 
quark models, symmetries and so on, can be applied to charmed hadrons. On the other 
hand, the presence of c-quark with mass of 1.5 — 1.8 GeV inside charmed hadrons allows 
one to use QCD methods which take this large parameter into account, like, for example, 
Effective Heavy Quark Theory (EHQT). 

At the moment the spectra of charmed hadrons, both mesons, and baryons are known 
much worse than the spectra of light hadrons. One of the reasons of this situation is due to 
a small lifetime of stable charmed particles. The development of new experimental methods 
based on the use of vertex detectors gives a hope to have a significant progress in study of 
charmed particle spectrum in the nearest future. 

The small lifetime of charmed particles makes it reasonable to look for charmed particles 
in experiments with their production in beams of high energy hadrons, when the lifetime of 
charmed particles in laboratory frame increases a few times. One of the goals of the new 
generation experiments with S~ beams to be performed at CERN (WA89) and at Fermilab 
(E791) is the study of charmed meson and baryon spectra. 

There is a wide variety of theoretical models which are more or less satisfactorily applied 
to description of charmed particle hadroproduction. These models can be selected into three 
main groups: the models based on perturbative QCD and parton model 0-0, the models 
of Monte-Carlo type based on Lund string model || and their subsequent development 
in the form of event generator PITHIA |J, and the models which incorporate the Dual 
Topological Unitarization scheme |H|-|[L0||. Among the last group we will distinguish the 
model of analytical calculation of spectra, so called Quark-Gluon String Model (QGSM) 
developed in ITEP |J-0i| . 



The use of perturbative QCD based models for calculation of charmed particle spectra 
leads to systematic underestimation of these spectra in the fragmentation region (x > 0.5). 
Different variants of Monte-Carlo models which produce a reasonable description of A c 
spectra at intermediate energies fail when applied at higher energies (see, for instance, 

The QGSM was successfully used for description of many features of multiparticle pro- 
duction, including inclusive spectra of secondary particles in hadron-hadron collisions both 



for light |[L1|-||15| and charmed [16|-|[22|| hadrons. In these papers the production of only 
stable particles (n, K, p, p, A, A c . . .) was considered. Spectra of meson resonances (p, / 2 , 
K*) in the framework of QGSM were analised in papers [^3, but the impact of reso- 
nances into 7r meson spectra was not considered in these papers and the contributions of 
direct production and of resonance decay were not separated. 

The main purpose of present paper is to apply the QGSM to calculation of charmed 
particle spectra produced in the beams of negatively charged hyperons, taking into account 
the decays of S'-wave resonances, like 1~ mesons (D* and D*), l/2 + (E c and S' c ) and 3/2 + 
(£*, S* and Q*) hyperons. The model parameters we use are obtained from the description 
of existing data on charmed particle production. The contribution of resonances into stable 
particle spectra will be taken into account in correspondence with their partial decay widths 

HJ. 
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2 Model description 



In the next section we will reanalize in the framework of the QGSM the inclusive spectra 
of stable charmed particles (A c , S c , Q c , D, D s ) taking into account the contributions from 
decays of corresponding resonances. 

We will calculate only the spectra integrated over transverse momenta p± of produced 
particles, although there exist modifications of the QGSM describing p± dependence p2j, |26j . 
S-wave charmed resonances decay into stable charmed particles with emitting 7r-meson or 



7-quantum [25 1, and we describe the kinematics of this decay as in [ O . 

Under such assumptions the invariant cross section of production of hadron h has the 
form 

da h da hdir J da R ^, , , 

x—— = x— h / x R - — ®{x R )dx R . (1) 

dx dx J dxR 

x*_ 

Here, x da ^ r is the direct production cross section of hadron h, and xr^^ is the i?-resonance 
production cross section. Function $(xr) describes two-body decays of resonance R into 
hadron h. After integration over transverse momenta of both hadron h and resonance R the 
function <&(xr) has the form 

/ x Mr 1 

In eqs. (|l|) and fl2|) xr is the Feynman variable of resonance R 



, M R x , M R x „ / 2yJ<pj>+m 2 

x +- E *_ p *> x -- E * +p *> ^-V X+X ±' x ±- ^ 

m is the mass of produced hadron h, Mr is the mass of resonance, E* and p* are energy 
and 3-momentum of hadron h in the resonance rest frame, < p\ > is the average transverse 
momentum squared of hadron h. 

The inclusive spectra of hadron h in the framework of the QGSM has the form ||11||-|T5|| 



00 



da h 

X-I- = J2°n(s)<p h n (x) , (3) 

ax n=0 

where cr n (x) is the cross section of n-pomeron shower production and <fn(x) determines the 
contribution of diagram with n cut pomerons. 

The expressions for cr n (s) and corresponding parameter values for pp and irp collisions 
are given in [llj [15]. For S^p interaction author of paper [2(| used the same values as for 



pp. In the framework of additive quark model we calculate the relation between pomeron 
residues in T,~p(E~p) and pp collisions 

7s p = 0.927 PP ; 7= p = 0.847 pp 
The function <fn{x) (n > 1) for irp interaction can be written in the form [PH]— 1|15|| 
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<p h n (x) = /*(*+, n)/*(x_,n) + /*(x+,n)/*(x_,n) + 2(n - l)/ a fc (x + , n)/ s h (x_, n) (4) 
and for baryon-proton interaction 



/nW = fl{x+,n)tf{x^n) + ft(x + ,n)fZ q (x_,n) + 2(n - l)/^(s + ,n)/£ a 0n_ ) n) , (5) 



where x± = | [w x 2 + x j± x] 



2 1 

eft/ 



The functions f^(x,n)(i = q, q, qq, q se a) hi @ and (H) describe the contributions of the 
valence/sea quarks, antiquarks and diquarks, respectively. They represent a convolution of 
quark/diquark momentum distribution functions Ui(x, n) in the colliding hadrons and of the 
function of quark/diquark fragmentation into hadron "h" G^(x,n): 

fi(x,n)=[ Ui(xi,n)Gi(x/xi)dxi . (6) 



The projectile (target) contribution depends on the valiable x + 

The quark functions f q (x,n) in the case of np and pp-collisions are given in |ll||-JiT 
without taking into account the possible resonance production. 

For the X - (H - ) beams the functions f^(x,n) are expressed in terms of corresponding 
s- (fg(x,n)) and d- (f$(x,n)) quark functions in the following form 

f^-\x,n) = \f^-\x,n) + §/f^(x,n) 

(7) 

f^-\x,n) = lf^~Kx,n) + lft {S ~\x,n) . 

In the framework of the additive quark model diquarks in S"-wave baryons may have 
spin (isospin) and 1. So the diquark functions fq q (x) are expressed in terms of scalar (0) 
and vector (1) diquark functions with the weights determined by SU (6) symmetric functions 



= + |yjg (x,n) + §/XM 

fT ] = \fT ] M + §/Jg>,n) + lf^\x,n) (8) 

f h q t ] = \f^~ ] M + §/gg>,n) + lf^\x,n) . 

In what follows, we will assume that distribution functions of scalar and vector diquarks 
u qq (x, n) are the same. Of course, different diquarks fragment into baryons in different ways. 
So, for instance, the direct production of A c in pp collision is determined by scalar (and 
isoscalar) diquark function f( u d) (Fig-l a ) ; an d direct production of S c and S* hyperons are 
determined by the vector diquark function f( u d)i- 

We also assumed that spin of diquark do not influence splitting of diquark (one of such 
terms are shown on Fig. lb). In the case of diquark fragmentation into meson this assumption 
leads to equation fM q ^(x,n) = fM q s{x : n) and eq.(H) reduces to the expressions 
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/£ W (*,») = \fuu {p) M + §^ (p) (*,n) 

f™V-){x,n) = \fT~\x,n) + |/^-)(x,n) (9) 



which coincide with p0| for the diquark in proton and £~-hyperon. 

A full list of the quark/diquark distribution functions in 7r-meson, p, E~, and H _ - 
hyperons used in this work is given in Appendix I. 

Further, we will assume that fragmentation functions of the quarks and diquarks does 
not depend on the spin of the picked up quark (or diquark). From this assumption the 
equality of the fragmentation functions of the corresponding quarks or diquarks to S c - and 
£*— , Sp- and H*-baryons, D- and D*-mesons follows. The method developed in the present 
paper for parametrization of the fragmentation function slightly differs from the one used 
in ||ll||-|pl||. We represent the fragmentation function as a sum of two terms. The first one 



is parametrized as a product of two polynomials each of which corresponds to the sum of 
all possible assymptotic and preasymptotic terms in the fragmentation (x — > 1, expanding 
in the series of (1 — x)) and central (x — > 0, expanding in the series of x) regions. The 
second term stands for the case when the fragmented object neither as a whole nor even 
partly gets into created hadron. In the case of nonleading fragmentation the corresponding 
functions are parametrized only using the second term. Appendix II contains the full list of 
fragmentation functions of quarks and diquarks into A c , £ c , £~, E c , E' c , S*, fl c , f2*-baryons 
and D, D*, D s , Z}*-mesons. The values of free parameters, determined from the comparison 



with the experimental data on A c production in np [29] and pp |30, |31] collitions, D and 



Z}*-mesons in pp |33|-||35|| and up p5|-p8| interation are also given in Appendix II 



3 Comparison with experimental data and the predic- 
tions of the model 

3.1 Description of charmed hadron spectra 

In this section, we consider the description of the existing experimental data for the 
A c -hyperon and D and D*-meson production in irp and pp-collisions in the framework of 
the present model. 

The A c -baryon spectra in 7r~p-collision at 230GeV/c are shown in Fig.2a,b corre- 
spondingly. The theoretical curves are calculated by using eq.(|I|) taking into account the 
contributions of A c produced in decays S c — > A c tt and S* — > A c ir. The dotted line in Fig. 2b 
shows the contribution of direct A c . As one can see the agreement with ir~p data is satisfac- 



tory. Concerning the data |3CJ, |3]J] in pp-collision one can see a noticeable difference between 
the data of groups |5D| and pi[ . This uncertainty does not allow us to have the unambigu- 
ous values of model parameters. The values given in Appendix II correspond to the curve 
plotted in Fig. 2. Our calculations of A c integral cross sections in irp and pp interactions are 
compared with corresponding experimental data in Table I. Large uncertainity in A c integral 
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cross section presented in paper [ 30 1 are mainly due to method used for extrapolation of the 
data to low x region. 

are 



The inclusive xp distribution of D-mesons produced in pp and irp collision p3|-p8 



compared with our calculations in Figs. 3-5. Both theoretical predictions for the integral 
cross sections and experimental data are shown in Table 2. 

In Figs.3a-d we plotted the experimental points for the spectra of different .D-mesons 
produced in pp-collision at 400 GeV/c j32] together with QGSM calculations. 



Inclusive distributions of all .D-mesons in pp interaction at momenta of 200 GeV [35 



400 GeV/c |32j and 800 GeV/c El] are presented in Fig.4a-c. The theoretical curves are 
calculated for the sums of spectra of all D-mesons. 

Fig.5a-d present a comparison of the x c -distribution of leading (D~ and D°) and non- 
leading (D + and D°) charmed mesons in n~p interaction at 200 GeV jc |35[ and 360 GeV '/c 



[37 1 with theoretical calculations. 

As far as we consider here for the first time the QGSM spectra of resonances taking 
into account their subsequent decays, it is important to compare our calculations with the 
available data on D*-meson productions in pp [3~2~| and irp collitions. The data on 
reactions ir~p — > D* + /D*~X and n~p — > D*°/D*°X at 360 GeV/c [|38| are compared with 
our predictions in Fig.6a,b. It seems that the agreement for the sum of spectra of D* + 
and D*~ mesons is rather reasonable. As to the neutral mesons D*° and D*°, there the 
experimental information is rather scarce. The curves in Fig. 6 correspond to the sums of 
spectra of D* + and D*~ (Fig. 6a) and D*° and D*° (Fig. 6b) mesons. Table III contains the 
experimental data on integral cross sections on D* meson production in pp and irp collision 
together with our calculations. The upper limit on .D^-meson integral cross section in pp- 
interaction 3^] is also shown in Table III together with our calculations. As one can see, 
agreement of model calculations for above mentioned inclusive spectra is rather good. 



3.2 Model predictions for charmed hyperon production with hy- 
peron beam 

In this section we consider the predictions for different inclusive spectra of charmed 
hadrons on S~ and S~ beams. We restrict ourselves to initial momenta of 340 and 600 GeV jc 
which correspond to the existing experiments in CERN (WA89) and FNAL (E791). The 
parameters used were obtained from description of experimental data on charmed hadrons 
in pp and 7rp collisions. The inclusive spectra of stable baryon states A c , H+, S° and fl c 
at 340 GeV/c are shown in Fig. 7a and the spectra of S* + and Q*° are plotted in 

Fig. 7b. The same spectra for 600 GeV/c are given in Fig. 7c, d. 

Fig.8a-d present our calculations for the same particles (except A c ) obtained in E~p 
collision. The predictions for Feynman-x distributions of spectra of D and D* mesons are 
shown in Fig.9a,b (340 GeV/c) and Fig.9c,d (600 GeV/c). The same calculations for E~ 
beam are presented in Fig.lOa-d. 

It is interesting to consider the predictions for D s and D* meson production in hyperon 
beams. As far as D~ and D*~ contain valence s-quark from initial hyperon beam, the 
cross-section of these mesons should not be small. Our calculations for the same momenta 
are pictured in Fig.lla-d (£~-beam) and Fig.l2a-d (S~-beam). 
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4 Conclusion 



In our paper a modification of QGSM which takes into account the resonance decay 
contributions is presented. The main result of our analysis of experimental data can be 
summarized as following: 

- The model under consideration describes with a good accuracy the inclusive x-spectra 
of charmed particles in np- and pp-collisions. So we hope that our predictions for charmed 
particle production of hyperon beams are fairly reliable. 

- The model predicts fairly large yield of resonances and as a result values of the 
spectra of unfavored stable particles in the fragmentation region are rather close to the 
spectra of favored particles. As a striking example we note the — > H+X reaction where 
at x F > 0.2 the spectra are completely determined by leading produced S*°-resonances 
via S*° — > 7r~ decay. 

Let us note that it is possible to generalise the present model and include the dependence 
on transverse momenta of produced particle taking into account the intrinsic charm and also 
to consider processes of production of P-wave resonances. In our opinion we have well 
described the data in the fragmentation region (x F > 0.2). More accurate description of the 
experimental data requires a more careful consideration of the resonance decay kinematics, 
that is essentially connected with taking into account the transverse momentum dependence. 

Acknowlegments. We are grateful to A.B.Kaidalov, K.G.Boreskov, M.A.Kubantsev, 
O.I.Piskounova, and for useful discussion. G.H.Arakelyan would like to thank the Theo- 
retical Department of ITEP for kind hospitality during his visits to ITEP where this work 
was done. This work was supported in part by International Science Foundation Grant N 
RYE000 and Grant INTAS-93-79. 
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Appendix I. The distribution functions of quarks and 
diquarks in the projectile and target hadrons 

The distribution functions of quarks (diquarks) in hadron h are parametrized in the 
standard form 

f?(x,n) = C i x ai (l-xfi , (A.l) 

where f3[ = + 2(n — 1)(1 — a®). The coefficients Cj in ( |A.l ) are determined by the 
normalization condition f™(x, n)dx = 1 and are equal to 

_ r(i + q)r(i + ^) 

Cj " T(2 + a + (3) ■ (A - 2) 

Here, r(a) is the Gamma-function. 

The values of a and (3 can be expressed in terms of Regge trajectory intercepts and are 
shown in Table A.l. We use the values of a Q p = 0.5, a° = 0, a% = —0.5. 

Appendix II. The functions of quark and diquark frag- 
mentation into charmed hadrons 

The fragmentation functions for both quarks and diquarks consist of two terms. The 
first one appears when the object i under consideration (quark, diquark as a whole or some 
its part) enters into the produced hadron h (favored fragmentation, see Fig. la). The second 
term corresponds to the absence of object % in the hadron h (unfavored fragmentation). 

G\ (x) = G h u (x) + Gi{x) (A.3) 
Functions G^x) and G%i(x) in the x — > 1 region contain the universal factor 

F^x) = (1 - x) x ~< (A.4) 

For convenience we will further consider the fragmentation functions of quark and di- 
quarks separately. 

II. 1 Quark fragmentation functions 

In the case of leading fragmentation functions G^x) are parametrized in the form 

G h u (x) = dV l F!(x)(l - xy k (A.5) 

where i = q(u,d),s is the type of the fragmented quark, k is the type of object (quark, 
antiquark or diquark) which combines with quark i to produce hadron h. The expressions 
for Si (lines 1, 2 for the case of quark fragmentation) and 7^ are given in Tables A.2 fl and 
A.3. Let us note that 7^ does not depend on the type of fragmented quark. 

1 Table A.2 contains also the values for diquark fragmentation. 
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The second term in (|A.3|) is parametrized in the form 



G* ( z ) = a%Fi(x)(l - x) Sm . (A.6) 

Here, m depends on the quark content of produced hadron h and the type of fragmented 
quark. 

In the case of nonleading fragmentation function Gf(x) contains only the second term 
given by eq.( |A.6j ). The 5 m is expressed in terms of intersepts of Regge trajectories as pre- 
sented in Table A. 4 (lines 1-3). 

II. 2 Diquark fragmentation functions 

Since only one quark from diquark takes part in meson creation, diquark fragmentation 
function into mesons is described by the same equations as the quark fragmentation function 
into baryon: eqs. (|A.3HA.6|) with the same values of 7^ and 5 m . In this case I — q,s are 
the quarks from diquarks which take part into meson creation, k = qjjc, sjjc (j denotes 
string junction) is the part of diquark width does not enter meson h and c-quark with string 
junction (see lines 3 and 4 in Table A. 3). 5 m , where m stands for sum of all components 
of fragmented and created objects, is given in lines 1-3 of Table A. 4. If only one of the 
quarks from diquark enters into produced meson, for instance G*£ (x) , G® s s (x) , (x) and 
G%t(x), eq. ([PD is divided by factor 2. 

As to diquark fragmentation into baryons, the functions G^(x) contain both leading 
term, described by diagram in Fig. la and terms, corresponding to the splitting of diquark 
(one of them is shown in Fig. lb). In this case G^x) is parametrized as 

Gf f (s) = (E^ £ " Oi + E ^i 1 - x^F^x) . (A.7) 

l k>l 

If fragmentation functions do not contain the leading term, G^x) has the form 

Gii(z) = (E^ £i ) • (E^ 1 -xV k )Fi{x) . (A.8) 

l k>l 

In eqs.( |A.7| ) and ( |A.8|) / means the type of objects from diquark % which take part in 
creation of baryon h (except the whole diquark) in x — > region. The values of E\ are taken 
from Table A. 2. k in ( |A.7|) and ( |A.8|) is the type of diquark splitting diagram (one of them 
is shown in Fig. lb) and denotes the remained part of diquark which do not form baryon h 
and the objects combined baryon h together with objects I. Table A. 5 contains the values 
of 7^ expressed in terms of of Regge trajectory intercepts for charmed baryon production. 

If only one quark or quark with string junction from fragmented diquark i passed into 
the created baryon, the corresponding constants bi and are divided by factor 2. 

The function G^x) for diquark fragmentation into baryon is also parametrized in the 
form ( |A.6|) . The corresponding values of 5 m are given in Table A.6. Here m is the sum over 
all types of objects in fragmented diquark i and created baryon h. The parameter values 
which were taken from the comparison with experimental data are equal to 



9 



a° = a°* = ao c = af c = a = 0.023 

Gq c = Oq c = CLq c = CLq " = CLq s = CLqO 

a\ = 4, a>2 = a 3 = 0,4 = d = 0.23, a 5 = . . . = a 13 = a 

61 = b 2 = 1, 63 = h = h = 0, h = b 7 = 0.1, 6 8 = 0.01 , 

5 = | stands for s-pair supression. Intercepts a° p) , a% are given in Appendix I. According 
to paper |39| we used only a° = —2.18 
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Table A.l 
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Table A. 2 
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Table A. 3 
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Table A.4 


N 


m 




1 




qqc 




2(1 -a p ) 




2 




qsc 




2 Q^p 




3 




ssc 




2(1 -<) 




4 




3qcj 




2(1 -a N ) 




5 


2qscj 


2(1 


— cxm) + ol p — 




6 


q2scj 


2(1- 


- a N ) + 2(a p - 


- Oiy) 


7 




3scj 


2(1- 


- a N ) + 3(a p - 


- a v ) 



12 



Table A. 5 
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Table A. 6 
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Table I 



Reaction 


Ref. 


P L (GW/c) 
or 




Vtheor&b) 


pp -> A c X 
all x 




30 




63 GeV 


40 ± 18 
204 ± 11 
2046 ± 836 


660 


pp^A c X 
\x\ > 0.5 




31 




63 Ge^ 


101 ±18 ±26 


84 


tt-N -> A c X 
x c > 




29 


] 


230 GeV/c 


4.9 ± 1.4 ±0.7 


6.8 
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Table II 



Reaction 



Ref. 



P L (GeV/c) 



O'theor(lJ'b) 



pp^D + X 



pp 
pp 
pp 
pp 



D~ X 
D° X 
D° X 
D+/D- X 



pp -> D°/D° X 

pp -> D+/D- X 

pp -> D°/D° X 

pN -> Z?/Z) X 
if > 



7T 



7T-X -> D+/L>° X 

n-p^ D+/D- X 
n-p -> X 



[0 




H 

a 



3|] 



35 



36, 37 



36, 37 



400 
400 
400 
400 
800 

800 

800 

800 

200 

200 
200 
200 

200 

360 
360 



5.7 ±1.5 
6.2 ± 1.1 
10.5 ± 1.9 
7.9 ±1.5 
33 ±7 

26 ± 14 



22 



-4 



1.5 ±0.7 ±0.1 



+0 4 
17 -0.3 ±0.1 



_i_0 5 
3 - 3 _0 4 ±0 ' 3 

+0 4 
2 ' 3 -0.3 ±ai 



4-fl 5 

3-2_ Q 4 ±0.2 



5.7 ±1.5 



10.1 ±2.2 
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Table III 



Reaction 


Ref. 


Pl (GeV/c) 




<?theor{nb) 


pp -> D* + /D*- X 


@ 


400 


9.2 ±2.4 


7.14 


pp^ D*°/D*° X 






400 


5.8 ±2.7 


8.8 


n-p^ D*+/D*- X 


[38 




360 


5 0+ 2 - 3 
-1.8 


5.0 


n-p^ D*°/D*° X 


[38 




360 


7.3 ±2.9 


4.5 


n-N -> D* + /D*- X 


[35 




200 


2.4 ±0.4 ±0.2 


2.6 


pp -> D+/D7 X 
if > 


[32 




400 


< 2.5 


2.8 



16 



Table captions 

I comparison of the experimental cross sections of A c production in pp and up inter- 
actions with results of our calculations. 

II The same as in Table I of D-mesons production. 

III Experimental data and model calculations for D*- and Z) s -mesons production in 
Tip and pp collisions. 



17 



Figure captions 

Fig.l Diagrams corresponding to fragmentation of diquarks into baryons. 

Fig.2 Comparison QGSM calculations with experimental data on A c spectra in: a) 
n-p -> A C X at Pi = 230 GeV/c f2|, b) pp -> A C X at ^ = 63 GeV/c 0, |TJ. 

Fig.3 The a^-distributions of different D-mesons in collisions at 400 GeV/c [Q: a) 
£>+ b) D~, c) Z?°, d) D°. 

Fig.4 Inclusive distributions of all D-mesons in pp interaction at a) 200 GeV/c |35|| , b) 
400 GeV/c H, c) 800 Ge^/c ||. 

Fig.5 comparison of the model calculations with experimental data on leading (D~/D°) 
and nonleading (D + /D°) charmed mesons in n~p interaction: a) leading, 
200 GeV/c H, b) nonleading, 200 GeV/c ||, c) leading, 360 GeV/c |7|, d) 
nonleading, 360 GeV/c 0- 

Fig.6 The x^-dependence of /}*-meson in 7r~p interaction at 360 GeV/c [ 3~8"| : a) 
D* + /D*- -mesons, b) L>* /.D* -mesons. 

Fig.7 Inclusive spectra of charmed baryons in S^p collision: a) A c , H+, 5°, f2 c at 
340 GeV/c; b) 5*°, 5*+, ft* at 340 GeV/c; c) the same as in a) but for 600 GeV/c; 
d) the same as in b) but for 600 GeV/c. 

Fig. 8 The same calculations as in Fig.7 for S~ beam. 

Fig.9 Predictions for different D and D*-meson production on X - beam: a) D + , D~, 
D°, D° at 340 GeV/c; b) D*+, D*~, D*°, D*° at 340 GeV/c; c) the same as in a) 
at 600 GeV/c; d) the same as in b) at 600 GeV/c. 

Fig. 10 The same as in Fig.9 for 5~ beam. 

Fig. 11 The x^-dependence of D s - and D*-mesons production in S~p collisions: a) Df 
and D~ at 340 GeV/c, b) D*+ and D*~ at 340 GeV/c, c) as in a) at 600 GeV/c, 
d) as in b) at 600 GeV/c. 

Fig. 12 The same as in Fig. 11 for S~ beam. 
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